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Abstract: We repeat and correct the recent calculation of the thermodynamic potential 
of hot QCD in the limit of large number Nf of fermions. The new result for the thermal 
pressure turns out to agree significantly better with results obtained from perturbation 
theory at small coupling. For large coupling, a nonmonotonic behaviour is reproduced, but 
the pressure of the strongly coupled theory does not exceed the free pressure as long as the 
Landau pole ambiguity remains negligible numerically. 
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1. Introduction 



In a recent paper p] the thermal pressure in QCD with a large number of fermions Nf 3> 
N c ~ 1 was calculated at next-to-leading order (NLO) in a large Nf expansion. Although 
the large-iVy limit is afflicted by the presence of a Landau pole, thermal effects can be 
studied in a cutoff theory provided the temperature is much smaller than the cutoff which 
in turn has to be smaller than the scale of the Landau pole. Then at NLO order of the 
large Nf expansion exact results, nonperturbative in the effective coupling g 2 s = g 2 Nf/2, 
can be obtained as long as g 2 s <C 6ir 2 . 

Exact large-TV results in scalar field theory at finite temperature have been obtained 
previously and used to study the (poor) convergence properties of thermal perturbation 
theory ||. An exact nonperturbative result for a more QCD-like theory is of particular 
interest in view of the various recent attempts to improve thermal perturbation theory 
in hot QCD by selective resummations [||, [f| g, 0, § , for which it may serve as a testing 
ground. In Ref. PD], it was proposed to interpret a failure of some technique at large Nf and 
reasonably large g 2 s as meaning that the technique is certainly not valid in full, small- Nf 
QCD. However, Peshier || recently argued that the strong-coupling behaviour of large Nf 
QCD is probably too different from that of small- Nf QCD to draw such conclusions. 

The result presented in Ref. [Q] is in fact very different from an ideal quasiparticle 
picture as pursued in Refs. [10, 11, 12 j. 1 According to Ref. Q, the gluonic contribution 
decreases as a function of g 2 s only up to a certain value of g 2 s , after which it rises and 
even exceeds the free pressure long before the coupling is so strong that the presence of a 
Landau pole becomes relevant. 

In the following, we shall present the numerical result that two of us (A.I. and A.R.) 
have obtained by a new implementation which closely follows the approach of Ref. yj . This 
result differs from that published by one of us (G.D.M.) in Ref. but after correcting the 
error in the computer code 2 underlying the latter, the two independent evaluations agree 
to an accuracy better than 1 : 10 . 

The new result turns out to follow rather closely the perturbative results to order g 5 
up to g 2 s ~ 6. At g 2 s ~ 12 the pressure goes through a minimum after which it rises, in 
qualitative accordance with the result presented in Ref. ||], but the exact result starts to 
exceed the free-gluon pressure only at values of g 2 s > 28, which is so large that the Landau 
pole starts to influence the results noticeably. 



2. Results 

The NLO contribution to the thermal pressure, of order Nf, is given by the one-loop 
gauge-boson contribution with any number of (renormalized) fermion bubble insertions 

1 A discussion of the HTL-quasiparticle picture of QCD thermodynamics underlying the approach of 
Ref. H in the context of large- Nf is contained in Ref. jl3| . 

2 In evaluating Eq. (3.11) of Ref. jjj], the imaginary part of the logarithm of the longitudinal propagator 
was calculated as the arctangent of the imaginary part over the real part without checking whether the 
argument was within the principal branch of the arctangent function. 
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Figure 1: Exact result for Pn lo / Pfree as a function of g^ffC/- 4 
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linear in g c s , in comparison with the previous result of Ref. |IJ] and two sets of perturbative results 
through order g 4 and g 5 : (a) with renormalization point chosen within a power of e of we fE T; 
(b) within a power of 2 of 2irT. The line marked "FAC" corresponds p, = ire^~~ lE T where the 
perturbative result to order g 4 coincides with the one to order g 5 . 



Q. Carrying out the sums over Matsubara frequencies, this expression involves terms 
proportional to the Bose distribution rib, which are best calculated in Minkowski space, 
and parts without this factor, which Ref. [IJ evaluated partly in Minkowski and partly 
in Euclidean space. To avoid spurious logarithmic divergences, it is crucial to employ 
a Euclidean invariant cutoff A when cutting out the Landau pole. This introduces an 
error which is suppressed, relative to the full thermal contribution, by ~ T 4 /A 4 , so that 
the ambiguity caused by the Landau singularity is well under control for A 3> T. If the 
coupling <C 67r 2 , the Landau pole is exponentially large and one may choose a large 
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cutoff A 2 = aAL andau , which following Q we shall vary by taking a between 1/4 and 1/2. 

To ensure Euclidean 0(4) invariance when performing parts of the calculation in 
Minkowski and parts in Euclidean space, which have to be connected by great arcs, one 
needs the analytic continuation of the complete fermion one- loop self-energy to the complex 
energy plane. The relevant formulae are listed in the Appendix. 

Ref. jjj calculated pieces linear in in Minkowski space. Terms without were 
computed along a complex frequency contour which ran up the Minkowski axis to uJ m3iX < 
ALandauV^ f° r some a < 1, then along the great arc to Euclidean space, and back down to 
qo = y Qmax ~ Q 2 \ finally, a Euclidean integration of the rib free term was performed over 
4-spheres in Euclidean space up to Q 2 < Aj^^a. 

It is in fact simpler to calculate all pieces linear in in Minkowski space, and all 
terms without rib in Euclidean space. By actually calculating both ways, we have a rather 
non-trivial numerical check on the result. In our numerical implementation both ways 
turned out to agree within numerical errors of about 10 -5 . 

In Fig. |] we give our numerical result as a function of 

g ^ s (fl = vre~ 7B T). 3 The new 

result agrees well with the perturbative results to order g 5 up to g 2 s ~ 5, where the 
renormalization scheme dependence of the y 5 -result is still reasonably small (the previous 
result of Ref. ]l] showed significant deviations from the perturbative results already for 
g 2 s <J 2). If the perturbative result to order g 5 is optimized by fastest apparent convergence 
(FAC), which requires that the result to order g 4 coincides with the one to order g 5 and 
which amounts to p, = 7re2 _7B T, the agreement with perturbation theory is improved and 
extends to g 2 s as 7. 

For higher values of g 2 s the exact result flattens out and reaches a minimum at g 2 s ~ 
12. For still higher values the pressure rises but, contrary to the previous result of it 
does not exceed the free pressure for the range of coupling considered in ]|] . 

In Fig. ^ we consider even higher values of g 2 s and find that eventually the thermal 
pressure grows larger than the free pressure. This occurs at g 2 s > 28 where ALandau/I" < 34. 
While this still seems to be a reasonably large number, the numerical result starts to become 
sensitive to the cutoff just where the pressure approaches the free one. The four curves 
displayed in Fig. |2| show the result of varying the parameter a in the UV cutoff -v/aALandau 
in the Minkowski and Euclidean parts of the calculation (om and ag, resp.) from a = 1/4 
to a = 1/2. The numerical result is rather insensitive to this below g 2 s ~ 25, but very 
sensitive in the region where the pressure starts to exceed the free one. 



3. Conclusion 



The exact result for the pressure of hot QCD in the limit of large Nf shows a nonmono- 
tonic behaviour as a function of the coupling. The minimum of the pressure is reached 
when ALandau — 480T and where the ambiguity introduced by the Landau singularity is 



Tabulated results can be obtained on-line from tittp : //hep . itp . tuwien. ac . at/~ipp/data/| 



In contrast to the exact result, the perturbative results depend on the value of the renormalization 
point ft, which we vary between ttT and AtyT, expressing everything as a function of g^sifi = 7re 1E T), 
however, to make a comparison possible. 
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Figure 2: The result for Pnlo/P{icc up to g^ ff = 35 and with different cutoffs. 



completely negligible. For higher values of the coupling the pressure eventually reaches 
and exceeds that of the free theory, but at that point the Landau pole is at ALandau < 34T. 
Above this point the result becomes increasingly sensitive to the precise cutoff which has 
to be chosen between T and AL an< j a u- This suggests that only the nonmonotonic behaviour 
is to be taken seriously, but not the fact that the free theory value is eventually reached 
and exceeded. 

So in contrast to the previous result of Q, the corrected one does not imply that 
an ideal quasiparticle picture (where the pressure has to be smaller than the free one) is 
necessarily in conflict with the actual physics of QCD in the limit of large fermion number. 
In order to be compatible with the nonmonotonic behaviour at large coupling, however, 
a quasiparticle picture would require a correspondingly nonmonotonic behaviour of the 
quasiparticle masses. While this is not particularly natural for the simple quasiparticle 
picture underlying the approaches of JTol, 11, O, this is not a priori excluded for the more 



complicated HTL-based ones of Ref. ]6[. This issue is discussed in more detail in Ref. [13]. 



A. Appendix: Gauge-boson self-energy 
A.l Spectral representation 

A convenient starting point for performing the analytic continuation of the self-energy from 
Minkowski to Euclidean space or vice versa is its spectral representation 

dpo 7i>0o,p) 



lWz,p) 



2tt 



Po 



(A.l) 



which is valid for any complex z [14|. The spectral form KuvijPOiP) is a purely real quantity 
that can be read off from the fermion loop evaluated in the imaginary time formalism 
according to 

d 3 k r°° dk r°° d qo 

(27)3 ^ 2^ 



(A.2) 
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xp (ko, k)p (go, q)-f ' ^ (A, q) 

with i 7 ^ = k^q u + g^/c 1 ' — g^ u k a q a + g^ u m 2 and q = k — p. The fermionic distribution 
function is given by rif(ko) = l/(e fco//T + 1) and the free spectral function by po(A;o,k) = 
27re(/c )(5(^-k 2 -m 2 ) = f^{d{k -ek)-d(ko+e k )) with e(k ) = k /\k \ and e k = Vk 2 + m 2 . 
(In the following we shall however consider only the ultrarelativistic limit m/T — > 0.) 

We need separately the transverse and longitudinal projection of the self-energy. Fol- 
lowing Weldon Q we define 2g 2 ff G = g^Tl^ and 2g^ s H = u^u^Yi^ with the thermal 
rest frame velocity = (1,0,0,0). 

Treating the various projections of the spectral density separately, we obtain the fol- 
lowing useful representations by analytically performing three of the four integrations: 

7Tx(f>0,P) = ^x(P0:P) ~ Kx(-P0,P) (A.3) 



„2 AT POO / 1 \ 

4(P0,P) = -^J o dk ("/(*) ~ 2J Ix (A.4) 
xe(k — po)0(\k — p\ < I k — p \ < \k + p\) 

where the ^-function stems from the angular integration between p and k (its usage here 
means #(true expression) = 1 and 6 (false expression) = 0) and X = G or H as in 

h = <rv(fc = *o =pI-p 2 (a.5) 

I h = u»u v I^{k = k) = ^(2k + p-p ){2k-p-p ). (A.6) 

The spectral density ir is manifestly real and odd in ko, i.e. 7r(fco,k) = — 7r(— A: ,k). 

To subtract the vacuum part, one just has to replace (nj(k) — ~) by nj{k). We shall 
do so in the following explicit results, because the vacuum part requires regularization and 
renormalization, after which the (Euclidean) self energy simply reads 

n ^c. = -|J (^ p2 ~ p " pu ) ( ln ^ - \ ) • (A - 7) 

A. 2 Minkowski result 

For Minkowski space we use the Feynman prescription 5 n F (po>p) = n(po + ipo^^v) f° r 
which the self-energy can be separated into 

Ren F ( Po ,p) = [°° ^^v)-^—r (A.8) 

J -00 2vr po - p' 

Imfi F (po,p) = ~l;e(p )ir(p ,p). (A.9) 

with P denoting the principal value as in = — — iir5(x). 

Note that with our expressions one has to turn the Euclidean po into the lower half of the complex 
plane po — > —v^ + e to obtain the retarded self-energy. 
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Inserting (|A.3|) in the expressions ( |A.8[ ) and ( A, 9 ) we reproduce the real part of the 
self-energy as given in Weldon's paper pi 



g 2 N f f 00 

Ren G (p ,p) = -rrir I dk n f( k ) 
/7T Jo 



4A; + ^_^!iog 

2p 



2k + p + p 2k - po + p 



2k + po — p2k — po — p 



(A.10) 



and 



ReIItf(j>o,p) 



9 2 N f 
2vr 2 J 



dk nf(k) 
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V v 
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4p 
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4p 
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(A.ll) 



The imaginary part was not explicitly calculated by Weldon, but we can provide a 
completely analytical result where no integration is left to be performed. It is given by 



ImIIx(po,p) 



-2 e{P0) ^p- 



? sApo+p\ 



\Po ~ P\ 



.„A/IP0+P| 

+e(Po + P)F X { — 7, — 



)-e{po-p)F x { 



(A.12) 



with symmetric and antisymmetric functions F x = (F x + F x )/2 and F x = (F x — F x )/2 
that are defined as 



F±(x) = I n f (k)I G (± P o,P,k)dk = (p z -p a )Fi( 



(A.13) 



2 2 

F*(x) = I n f (k)I H (±po,p,k)dk = ^^F 1 {x)T2poF 2 {x) + 2F 3 (x), 



where the Fi{x) are the following integrals 

F x (x)= / n f (k)dk = -x + T\og(e x/T + 1), 

J X 



(A.14) 



Fo(x) 



F*(x) 



tt 2 T 2 r 2 

kn f (k)dk = —— - — + xT\og(e x ' T + 1) + T 2 U 2 (-e x / T ), (A.15) 
6 2 

k 2 n f (k)dk = ~ + x 2 T\og{e x/T + 1) 



+2xT 2 U 2 {-e x l T ) - 2T 3 Li 3 (-e a; / T ), 



(A.16) 



with Li n (x) being the polylogarithm function. Note that Fq = simplifies our expression 
for Indie considerably. 
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A. 3 Euclidean result 

For Euclidean space we set z = iui and (using the antisymmetry property of the spectral 
density) we get 

tEucW ■ , „% _ f°° d Po„, , -pa ~ Eucl 



Ren tucl (^,p) = / -^n(p Q ,p)^j^, ImlT^(^P) = 0. (A.17) 

2-7T UJ 1 + J?q 



oo 



We are left with real integrals of the form 

/ dp Q ^^I G {pv,p 1 k) = -p 2 + {u 2 +p 2 )\og{uj 2 +p 2 ) (A.18) 

J LO A +Pn 



and 



/ 



dpo 9 2P ° 2 i h(po,P, k) = ^Po(8k - p ) - AkuJ arctan (—) 
uj z +Pq 2 \uj J 



- l -{Ak 2 -p 2 -uj 2 )\og(uj 2 +p 2 ). (A.19) 



With the appropriate integration limits we finally obtain the self-energy in Euclidean space 
as 

^ rt-^jf - »,(*> (« + ^ I. ) (A,0) 

g 2 N f f°° 
ReIL H (iuj,p) = ^2 J dknj(k) 



- 2at - r00 [ to 2 +p 2 -4k 2 uj 2 + (2k-p) 2 
2k H log ■ 



Ap fa u 2 + (2k +p) 2 



2ku ( 2k -p p 2k+p\ 

I arctan h 2 arctan arctan 

p \ UJ uj uj ) 



(A.21) 



This is in principle the result given in the Appendix of Q where the three terms involving 
the arc tangents are replaced by a common logarithm according to 

. 2 k-p p 2k + p i 1 + (u-i V y , . > 

arctan h 2 arctan arctan = log -p> — . A. 22) 

uj uj uj 2 & i + T Jhl- y ' 

However, while taking the principal branch of the arctan functions gives a smooth function 
over all k, on the right-hand side one must not restrict to the principal branch of the 
logarithm. 

For verifying the path independence of the numerical results we also need the self 
energy for complex energies. These may be obtained either from the analytic continuation 
of the results ( |A.2C| ) and ( A.21[) or from the spectral representation according to 



n(o + ib, p) = f ^-ir(p ,p) (- — ° 2 - i-, \ 2 2 ) , (A.23) 

J 2vr \(a-p ) 2 + b 2 (a - p ) 2 + b 2 ) 

with real and unambiguous integrals (for 6^0). 
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